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Quantum-state engineering, i.e., active manipulation over 
the coherent dynamics of suitable quantum-mechanical sys- 
tems, has become a fascinating prospect of modern physics. 
Here we discuss the dynamics of two interacting electrons in 
a coupled quantum dot driven by external electric field. We 
show the two quantum dots can be used to prepare maximally 
entangled Bell state by varying the strength and duration of 
an oscillatory electric field. Different from suggestion given 
by Loss et ai[Phys. Rev. A, 57 (1998) 120], the present 
entanglement involves the spatial degree of freedom for the 
two electrons. We also find that the coherent tunneling sup- 
pression discussed by Grossmann et aZ.[Phys. Rev. Lett., 67 
(1991) 516] persists in the two-particle case, i.e., two electrons 
initially localized in one dot can remain dynamically localized, 
although the strong Coulomb repulsion prevents them behav- 
ing so. Surprisingly, the interaction enhances the degree of 
localization to a larger extent compared to non-interacting 
case. We can call this phenomenon Coulomb-enhanced dy- 
namical localization. 

PACS numbers: 73.23.Hk, 73.61.-r, 78.66.-w, 78.7.+p 
Keywords: coupled quantum dot, dynamical localization, en- 
tanglement 



I. INTRODUCTION 

Coherent control of quantum systems has been at- 
tracted considerable attentions in recent years. A basic 
ingredient of quantum control is field-induced localiza- 
tion of a single electron in a double-trap [j] . The initial 
efforts devoted to acquiring conditions to maintain ex- 
isting localization with oscillatory electric field, and to 
create and maintain localization with a semi-infinite os- 
cillatory field d|| ■ It was found that the investigation of 
localization can be approximated by a two-state model 
consisting of the lowest symmetric and antisymmetric 
states of the double-trap potential In this case, 

perfect localization can be achieved with a strong time- 
periodic electric field that causes the Floquet quasiener- 
gies to be degenerate. Later, localization in supcrlatticc 
systems Q, in dissipative environments |7|, in molecu- 
lar systems induced by ultrashort laser pulses, and 
in trapped Bose-Einstein condensates HJl(| > by means °f 
oscillatory magnetic fields, has been studied. 

When two or more interacting particles are present, 
apart from the highly nontrivial problem of whether the 
strong many-body interaction can be overcome for the 
particles to create and preserve localization, the possi- 
bility of entanglement of the many-body wave functions 



arises |11|,|12[. Entanglement is an essential ingredient in 
any scheme of quantum information processing like quan- 
tum information cryptography and quantum computa- 
tion, and therefore it is a problem of great current inter- 
est to find or design systems where entanglement can be 
manipulated Jl3| . Most of the theoretical and experimen- 
tal activity until now has been associated with atomic 
and quantum-optic systems. Two- |14 1j|, three- [fj6| , 
and four-particle p7[ entanglement have been success- 
fully demonstrated experimentally in trapped ions, Ryd- 
berg atoms, and cavity QED. However, a further increase 
of the number of entangled particles in these systems is 
expected to be a severe experimental challenge. Recently, 
solid-state realizations of the entanglement have received 
increasingly attention due to the fact that semiconductor 
nanostructures such as quantum dots (QDs) and double 
quantum dots (DQDs) with well-defined atom-like and 
molecule-like properties, have been fabricated and stud- 
ied by many groups jl8|,[l9]]. Kane has proposed a 
scheme which encodes information onto the nuclear spins 
of donor atoms in doped silicon electronic devices where 
externally applied electric fields are used to perform log- 
ical operations on individual spins. Loss and DiVincenzo 
pl[ have presented a scheme based on spin exchange in- 
teraction effects. More recently, Imamoglu et al. j2jj have 
considered a quantum computer model based on both 
electron spins and cavity QED which is capable of realiz- 
ing controlled interactions between two distant quantum 
dots. Quiroga and Johnson |24|] have suggested that the 
resonant transfer interaction between spatially separated 
excitons in quantum dots can be exploited to produce 
maximally entangled Bell states. 

In the present work we study the coherent control of 
the quantum system consisting of two interacting elec- 
trons in a coupled quantum dot (see Fig. 1). With the 
initial state chosen to be in the spin-singlet space, the dy- 
namics is reduced to be confined to a three-dimensional 
Hilbert space, in which the three basis vectors are equiva- 
lent to the eigenstates of z component of spin-1 operator. 
We show that the maximally entangled Bell state can 
be prepared and maintained with a pulse of oscillatory 
electric field. Also we find that although the Coulomb 
repulsion between the two electrons is very strong, dy- 
namical localization can fully build up in the system pa- 
rameter manifold which corresponds to the exact cross- 
ing of the quasienergies developed from the unperturbed 
nearly-degenerate levels. 
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II. THE MODEL 

The Hamiltonian which we use to describe the dynam- 
ics of two interacting electrons in a coupled quantum dot 
driven by electric fields is 



H(t) = h(n, Pi ,t)+C, 



(la) 



i=l,2 



h(r, p ,t) = ^-- ezF(t) + Vi (r) + V v (r) , (lb) 
2m 



C 



K\r 1 - r 2 



(lc) 



where C is the Coulomb interaction and h the single- 
particle Hamiltonian. The dielectric constant k and the 
effective mass m are material parameters. The poten- 
tial Vi in h describes the lateral confinement, whereas 
V v models the vertical double-well structure. The lateral 
coupling of the dots is modeled by a quartic potential 



Vi(x,y) 



: 2/ 2 i 2\ 

-a (x + y ), 



(2) 



where we have introduced the isotropy parameter a de- 
termining the strength of the vertical relative to the 
lateral confinement. The lateral effective Bohr radius 
a-B// = y/fi/(muj z a) is a measure for the lateral exten- 
sion of the electron wave function in the dots. It has been 
shown in experiments with electrically gated quantum 
dots in a two-dimensional electron system that the elec- 
tronic spectrum is well described by a simple harmonic 
oscillator |^5) . In describing the confinement V v along the 
inter-dot axis, we have used a (locally harmonic) double 
well potential of the form 



V v — ■ 



2\2 



(3) 



which separates into two harmonic wells of frequency 
uj z (one for each dot) in the limit a ^> clb±, where a is 
half the distance between the dots and obi = y/ H/itiuj z 
is the vertical effective Bohr radius. Although in princi- 
ple a square-well potential would be a more accurate de- 
scription of the real potential than the harmonic double 
well, we believe there is no qualitative difference between 
the results presented below obtained with harmonic po- 
tential and the corresponding results using square-well 
potential. Finally, the time-dependent electric field has 
a DC- AC form, i.e., F(t) — F + Fi sinuit. 



III. SPIN-1 REPRESENTATION OF THE 
HAMILTONIAN 

Following Burkard et al.[22] we employ the Hund- 
Mulliken method of molecular orbits to describe the low- 
lying spectrum of our system. This approach accounts for 



double occupancies and is therefore suitable for investi- 
gating the questions at issue here. As a starting point for 
our calculations we consider the problem of an electron 
in a single quantum dot. The one-particle Hamiltonian 
by which we describe a single electron in the right (left) 
dot of the double-dot system is 



P 

2m 



+ ^( a V + y 2 ) + (z T a) 2 ), (4) 



and has the ground-state solution 



¥>±a(r) 



Tlfl 



3/4 



'aexp 



(5) 



corresponding to ground-state energy e± = hw z {l + 
2a) 1 2. The two local ground states are not orthogonal 
and their overlap is 



S 



d 3 np* +a (r)ip- a (r) = exp(-d 2 ), 



(6) 



where the dimensionless parameter d = a/as± denotes 
the distance between the two dots. A nonvanishing over- 
lap S implies that the electrons can tunnel between the 
dots. From these non-orthogonal states, we construct the 
orthonormalized one-particle wave functions 



$+ Q (r) 



*-.(r) 



1 



y/l - 2Sg - .g 2 



y/1 - 2Sg - g 2 



(cp +a - gip-a), 



(<P-a ~ gtP+a), 



(7a) 



(7b) 



with g = (1 — y/l — S 2 )/S. For appropriate values of 
system parameters such as the inter-dot distance, the 
overlap S becomes exponentially small as given in Eq. 
(7). In this limit an electron in one of the states $+ a , 
^-a is predominately localized around (x, y, ±a). In the 
following we consider this case and use these states to 
define the localized single-particle state. Schlicmann et 
al. J|^] have used theses two local states as qubits, which 
are realized by the spin state of an electron in cither 
orbital $+ a , or orbital $_o- 

Using $± a we generate six basis functions with respect 
to which we diagonalize the two-particle Hamiltonian 
H: The states with double occupation, ^^(ri, r 2 ) = 
$ T a(ri)$ Ta (r 2 ) [denoted by (1,0,0,0) T and (0,1,0,0) T , 
respectively] and the states with single occupation 
*i(ri,r 2 ) = [$ +0 (n)$_ a (r 2 ) ± *_ (n)$ +0 (r 2 )] [de- 
noted by (0,0,1,0) T and (0,0,0,1) T , respectively]. Calcu- 
lating the matrix elements of the Hamiltonian H(t) in 
this orthonormal basis we get 



H{t) 



H^t) 

{2e + U 3 )I lxl 



(2e + U 3 )I lxl ®H 1 (t), 
(8) 
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where 7i x i is a 1 x 1 unit matrix, and Hi is 

2e + Ui + fj,{t) V2w v 
ffi(i) = ( V2w 2e + U 2 V2w 

v y/2w 2e + Ux- n(t) 



(9) 



with the matrix elements given by 

e = (*± a \h(z T a)|* ±0 ), C/i - (*L|C|*i ), (10a) 



t/ 2 = (*;|C|*;),cr 3 = (* s _|c|* s _), 



(10b) 



™ = (*±a\h{z T «)l$ T a> + (n|C|*±a)^ = <*+a|C|*! >, 

(10c) 

. / n 2ea(l - g 2 ) „ . 

H{t) = Mo + Mi sin(wt) = - _ 2 - ^ (F + ^i sinwt). 

(lOd) 

Clearly, i/i denotes the intra-dot Coulomb interaction, 
whereas J/2 and C/3 describe the inter-dot Coulomb in- 
teraction, without particle transfer, w describes single- 
particle tunneling induced by dot-dot coupling and 
Coulomb interaction, pn = 2eFia(l — g 2 ) /(l — 2Sg + g 2 ) 
(i = 0, 1) describes the coupling strength between the 
electrons and external electric fields. 

Obviously the spin-triplet state is an eigenstate 
of the Hamiltonian (8) with the electron number on 
each quantum dot invariably one and has no response 
to the presence of the electric fields. This is one im- 
portant fact involving spin-triplet space. Hence we will 
focus our attention on the reduced spin-singlet Hamilto- 
nian H\(t). Furthermore, the matrix element v is the 
amplitude of simultaneous transfer of two electrons and 
therefore denotes cotunneling process. Compared to the 
single-electron tunneling term w, this term can be negli- 
gible. Thus the reduced Hamiltonian can be conveniently 
rewritten in terms of spin-1 operators 



Hi(t) = (2e + U 2 ) - n{t)J z + uJ 2 + 2wJ x 



(11) 



where u = U\ — U 2 is the effective Coulomb interaction, 
and Ji (i = x,y, z) are spin-1 operators defined as 



1 



J T, , — 



.,101 
v 2 \ 1 



0-10 

"7= I 1 _1 
v 2 \ 1 



(12a) 



(12b) 




(12c) 



and *f?+ a to the state 
whereas the delocalized state is identical 



Therefore, the localized two-particle state \&1 is equiva 
lent to the eigenstate | j z = 1) of J. 

Ii, = -i) 

to the state | j z =0). In the following we will denote 
with I LL) , with I RR), and V s + with | LR). 
The first term in Eq. (11) denotes a constant energy 
shift, and will be neglected in the following discussions. 
The evolution of any initial state | VE'(O)) under the action 
of H\ in Eq. (11) can be expressed as | *&(t)) = C\(t) \ 
LL) + C 2 (t) I LR) + C 3 (t) I RR) where the coefficients 
Ck (t) are determined by the time dependent Schrodingcr 
equation (In the following we set H — 1) 

u + fi(t) V2w \ / d 

y/2w V2w C 2 
V2w u - fj,(t) J V C * 

(13) 

and the chosen initial condition | ^(O)). In the absence 
of external electric fields, the eigenenergies of Hi can be 
easily solved as 




E 1 = ^(u- Vu 2 + 16w 2 ) 



E 2 = u, 



E 3 = i (u + yju 2 + 16w 2 ) , 
and the corresponding eigenstates are given as 



\v[ S) ) = ^{\LL)--^-\LR) + \RR)), 



W ( 2 A) ) 



— (\RR)-\LL)), 



^ 3 S) > = y{\LL) 



V2w 



\LR) + \RR)), 



(14a) 



(14b) 



(14c) 



(15a) 



(15b) 



(15c) 



and 



w here we ha ve defined the nor malization constants X = 
^jAw 2 + E 2 /V2w and Y = Aw 2 + E 2 /V2w. The su- 
perscript S (A) on the left sides of Eq. (15) denotes sym- 
metry (antisymmetry) under the spatial reflection. We 
can see from Eq. (15) that due to strong Coulomb re- 
pulsion, the symmetric ground state is dominated by the 
delocalized state \LR), whereas, the other two eigenstates 
are nearly degenerate and dominated by the two local- 
ized states \LL) and \RR). Note that although \ip 2 A ^) and 
l^g 5 ') look like a doublet in a single-electron double-trap 
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system which consists of a pair of symmetric and anti- 
symmetric single-particle states, there are fundamental 
differences for the present two-particle system. In fact, 
the superposition of the two localized states \LL) and 
\RR) implies that the spatial wave functions of the two 
electrons have been entangled, in the usual sense that 
they are not factorized into single-particle states. To de- 
scribe the degree of entanglement, we define the maxi- 
mally entangled Bell state 

\^ Bell ) = ^=(\RR)+e^\LL)), (16) 

with arbitrary phase angle <fi. Therefore the probability 
PBeii for finding the maximally entangled Bell state in a 
coupled quantum dot is given by 

PBeU = \\C 3 {t)+e i +C 1 (t)\ 2 . (17) 

In the following we will investigate the various energy 
spectrum and consequent dynamics for a GaAs (m = 
0.067m e , k = 13.1) system comprising two equal dots 
with vertical confinement energy uj z = 16meV (ag — 
36nm) and lateral confinement energy au> z = 8mcV. The 
inter-dot distance is chosen to be a = 20nm. The cor- 
responding system parameters are calculated to be u = 
5.6meV, and to = — 0.15mcV, implying \ip\ ) ~ \LR), 
and ~ (1/V2)(\RR) + \LL)). 

IV. LOCALIZATION PREPARATION FROM 
DELOCALIZED GROUND STATE 

We start the search for localization with the simplest 
case, a constant electric field F(t) = Fq. Before t = 
the system is in the delocalized ground state | tp[ S '), and 
at t = the field F is switched on suddenly. We find 
that when the value of electric field satisfies the following 
resonance condition 

F = u, (18) 

the initial delocalized ground state may be excited into 
a localized state with two electrons occupying the same 
right dot. To elucidate this resonance property, we show 
in Fig. 2 time evolution of the probabilities PLii(t) — 
|C 2 (t) | to find the two electrons in the different dots 
(solid line), Phh{t) = I Ci (*) 1 2 to find the two electrons in 
the left dot (dashed line), and Pr.r(£) = | C3 (t ) | 2 to find 
the electrons in the right dot (dotted line). It shows in 
Fig. 3 that both PLn(t) and -Pr_r(£) oscillate between 
and 1 with a definite period, while Pll^) always negli- 
gibly small during time evolution, meaning that a com- 
plete resonance takes place between the delocalized state 
I LR) and localized state | RR). In this case, because 
the population of the localized state | LL) in the left dot 



remains almost zero during time evolution, we can ne- 
glect its contribution and describe the dynamics by the 
reduced Schrodinger equation 

i ±(C 2 {t)\_(0 V2w\(C 2 (t)\ 
l dt{c 3 (t) )-{V2w J [C 3 (t) )■ 

Thus with the initial state | *(0)> =| ip[ S) ) ~| LR), we 
have the evolution of the system as follows 

C 2 (t) = cos(V2wt), (20a) 

C 3 (t) = cxp(in/2) sm{V2wt). (20b) 

Clearly our two-state approximation Eq. (20) describes 
the dynamics very well when compared with the exact nu- 
merical result shown in Fig. 2, implying complete Rabi 
oscillation between the localized state | RR) and delocal- 
ized state I LR) with oscillation period n / (\/2\w\) . 

Once the two electrons are localized, they can be forced 
to stay localized permanently by switching the field to 
another nonzero value [see Fig. 3(a)]. We show this effect 
in Fig. 3(b). This trivial way of keeping localization 
in a two-electron double-dot system originates from the 
fact that, due to the presence of constant electric field, 
the consequent energy mismatch among the three two- 
particle states prohibit tunneling from the localized state 
I RR) to the other states. 

Another way to implement localized two-particle state 
is through adiabatically varying the constant electric 
field, which induces a series of avoided crossings in the en- 
ergy spectrum. Consequently, the electron number in the 
ground state will experience a series of Coulomb stairs. 
This field-induced adiabatic localization is shown in Fig. 
4, where Fig. 4(a) plots the evolution of electron num- 
ber in the right dot for the ground state as a function of 
the constant electric field, and Fig. 4(b) plots the corre- 
sponding energy spectrum. It reveals in Fig. 4 that on 
adiabatically increasing the strength of constant electric 
field, a series of avoided crossings develop in the energy 
spectrum. Consequently, quantum transition occurs at 
these avoided crossings, resulting in a series of Coulomb 
stairs. We notice that adiabatically increasing the elec- 
tric field to a typical value of 2.8kV/cm enables complete 
localization of the ground state. 

V. ENTANGLEMENT OF TWO ELECTRONS IN 
THE PRESENCE OF CONSTANT ELECTRIC 
FIELD 

As shown in Fig. 2, a resonant constant electric field 
can induce a complete oscillation between the delocal- 
ized state \LR) and localized state \RR). If the constant 
electric field is turned off [see Fig. 5(a)] at time when 
the two electrons are fully localized in the right dot, as 
shown in Fig. 2(b), the strong Coulomb repulsion will 
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induce the resonance between the localized states | RR) 
and | LL) during subsequent evolution, whereas the de- 
localizcd state | LR) is inhibited to be occupied. This 
dark property of the delocalized state is shown in Fig. 
5(b). It reveals in Fig. 5(b) that the value of Plr is al- 
most zero during time evolution, suggesting that the two 
electrons are never separated into different dots. While 
cycling from one dot to the other, the two electrons are 
always correlated and entangled, and very likely to be 
found in the same dot. 

The entanglement between the two electrons illus- 
trated in Fig. 5(b) can be well described by a two- 
state approximation. Because the population of the 
delocalized state | LR) remains very small after time 
t = 7r/(2\/2|H) shown in Fig. 4(b), we can approxi- 
mate C2(t) (t > to) in Eq. (13) to first order of w/u 

C 2 (t) = exp(-i«t)[C 1 (t) + C 3 (t)}. (21) 

u 

By introducing C^it) from Eq. (21) in the Schrodingcr 
equation we reduce the system to an effective two-level 
system. The reduced equation has the form 

dt { c 3 (t) )-{ -2 W y u u )\ c 3 (t) J ■ w 

Thus with the initial state | *(t )) = -i \ RR) [see Eq. 
(20)], we have the following time evolution of the system 

d(t) = exp(-mt) sin{2w 2 t/u), (23a) 

C 3 (t) = -icxp(-iut)cos(2w 2 t/u). (23b) 

Substituting Eq. (23) into Eq. (17) we have the prob- 
ability to find the Bell state (| RR) + | LL))/\/2 at 
time t 

PBeii (t) = \ [1 + sin(ovt) cos(0 + tt/2)] , (24) 

where ui r = 4w 2 /k. In particular we can see from Eq. 
(23) that the system's quantum state at time 

r = itu/8w 2 +t Q , (25) 

corresponds to a <fi = —tt/2 maximally entangled Bell 
state (| RR) -i | LL))/V2. 

We present in Fig. 5(c) (solid line) the probability for 
finding the maximally entangled Bell state ((f) = —tt/2) as 
a function of time. The system parameters are the same 
as that used in Fig. 5(b). Clearly our two-state approx- 
imation [Eq. (24)] describes the system's evolution very 
well when compared with the exact numerical solution 
shown in Fig. 5(c), implying that the system's quantum 
state at time r corresponds to a maximally entangled 
Bell state (cf> = —tt/2). However, the degree of entan- 
glement degrades after time t, a consequence of the fact 



that the state at time t is not an eigenstate of the field- 
free Hamiltonian. So a single pulse of constant electric 
field can not preserve the entanglement in our system. 
Note that no maximally entangled Bell-state generation 
is possible if the effective Coulomb interaction k, along 
with the electric field, is turned off, as shown in Fig. 
5(c) (dotted line). This implies essential role of the non- 
linear Coulomb interaction in forming the entanglement 
between the electrons. 



VI. ENTANGLEMENT OF TWO ELECTRONS IN 
THE PRESENCE OF OSCILLATORY ELECTRIC 
FIELD 

We turn now to discussion of the entanglement in the 
presence of a sinusoidal field of the form F(t) — F\ smujt. 
Figure 6 illustrates the spectral features by plotting the 
Floquet quasienergies as a function of driving frequency 
u>. The amplitude value of the field is chosen to be 
F\ = 1.5kV/cm. It shows in Fig. 6 that on adiabati- 
cally increasing the value of w, the quasienergies e\ and 
e 3 approach each other. Especially when the value of 
driving frequency satisfies 

lo = u, (26) 

an avoided crossing occurs in the quasienergy spectrum. 

To elucidate the effect of the avoided crossing displayed 
in Fig. 6 on the quantum mechanical behavior of the 
system, we examine the dynamics of the system starting 
with unperturbed ground state, i.e., 1^(0)) = \ip[ S>> )- Fig- 
ure 7(a) shows time evolution of the probabilities PLR.{t) 
(solid line), PLL(t) (dashed line), and Pnjt(t) (dotted 
line), With the system parameters corresponding to the 
avoided crossing shown in Fig. 6. It reveals in Fig. 7(a) 
that the two electrons oscillate between the delocalized 
state | Li?) and two localized states \LL) and \RR). The 
occupations of two localized states are always the same 
and the oscillations are in-phasc. The maximum values 
of Phhit) and Pflfl(t) are 0.5, which corresponds to zero 
occupation of state \LR). Figure 7(b) shows the prob- 
ability pBeii to find the maximally entangled Bell state 
with (j) — tt. We can see from Fig. 7(b) that the degree 
of entanglement varies with time. In particular when 
the occupation of delocalized state \LR) is zero, the two 
electrons are maximally entangled with pseii = 1- 

Once the two electrons are in the maximally entan- 
gled Bell state, they can remain maximally entangled by 
suddenly turning off the oscillatory electric field [see Fig. 
8(a)]. We show this effect in Figs. 8(b)-(c) where Fig. 
8(b) plots time evolution of the occupations of three two- 
particle states and Fig. 8(c) the probability to find the 
maximally entangled Bell state with (f> = tt (solid line) . It 
shows in Fig. 8(b) that a pulse of oscillatory electric field 
induces the two electrons to stay on the same dot, while 
each of them occupies either of the dots with the same 
probability. It reveals in Fig. 8(c) (solid line) that the 
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two electrons remain maximally entangled during time 
evolution after the oscillatory field is turned off. This is 
different from what is shown in Fig. 5(c) where the degree 
of entanglement varies with time. Therefore, by a pulse 
of oscillatory electric field with frequency satisfying the 
resonant condition Eq. (26), the maximally entangled 
Bell state can be created and maintained. Loss and Di- 
vincenzo studied entanglement in double quantum dots 
involving the spin degree of freedom pj| . Here we have 
identified a complementary method that creates and pre- 
serves entanglement between the spatial wave functions 
of two electrons in a coupled quantum dot. 

The generation of maximally entangled Bell state 
shown in Fig. 8 can be well described by a two-state 
approximation. Taking into account symmetric proper- 
ties of the three unperturbed eigenstates of the system, 
the dynamics is determined by the one-photon transition 
between the ground state |<^ 5 ^) and the first excited two- 
particle state \<p^), whereas the transition from \(p\ S>> } 

to \(p 3 ) is prohibited due to their identical symmetry. 
In this case, we can approximate the Hamiltonian Hi(t) 

in a Hilbert space spanned by the states \ip± ^) and \tp 2 A ^) 



T = TTU/ {2wjJL\), 



(31) 



Hi(t) 



■V2m(t)/x 



-V2fii(t)/X 
E 2 



(27) 



In the interaction representation and after applying 
the rotating- wave approximation (RWA) the Schrodinger 
Equation takes the form 



_ d i 
" dt \ d 2 







a" 



(28) 



where d\ and d 2 are the probability amplitudes of the 
bare states |^>i ) and |<£>2 )j = V\/\/2X is Rabi fre- 
quency, and A = E 2 — E\ — uj is the detuning of the 
driving frequency uj from the transition frequency. In 
the case of one-photon resonance E 2 — E\ = uj, we ob- 
tain the time evolution of initial ground state [d\ (0) = 1 , 
^2(0) = 0) as follows 



di 



cos il r t, 



d 2 = i sinf2 r i. 



(29a) 



(29b) 



Thus the system oscillates between the ground state 
wherein two electrons are highly delocalized and the first 
excited state wherein the two electrons are highly entan- 
gled. Substituting Eq. (29) into Eq. (17), and in the 
weak coupling limit w <C u, we obtain the expression 
for the probability to find the maximally entangled Bell 
state 



pBell{t) = -(1 



1 sin 2 (w /lit /u), 



(30) 



where we have approximated Rabi frequency f2 r with 
w\x\ju. In particular we can see from Eq. (29) that 
the quantum state of the system at time 



corresponds to a <p = tt maximally entangled Bell state 
(\RR)-\LL))/y/2. 

The result of Eq. (30) is shown in Fig. 8(c) (dotted 
line). Clearly, in comparison with the exact numerical 
solution, our two-state approximation describes the sys- 
tem evolution very well, suggesting that the quantum 
state of system at time r = iru/(2w(ii) corresponds to a 
maximally entangled Bell state of the desired form with 
(f) = 7r. Therefore we arrive at the conclusion that a 
selective pulse of oscillatory electric field with duration 
t = iru/(2wfii) can be used to create and maintain a 
maximally entangled Bell state ((f> = it) in the system of 
two electrons in a coupled quantum dot. 

We notice from Eq. (31) that Bell-state generation 
time is significantly shortened by increasing the ampli- 
tude of the oscillatory electric field. This is important 
because shorter Bell-state generation time is fundamen- 
tal to the experimental observation of such maximally 
entangled state, which is impeded by inevitable decoher- 
ence occurred in the realistic double quantum dot sys- 
tem. The decoherence is the most problematic issue per- 
taining to most quantum computing processing. In the 
present entangled state proposal, the decohering time de- 
pends partly on the fluctuation of the single particle en- 
ergy caused by the modification of the confining poten- 
tial due to phononic excitations. There is also a quan- 
tum electrodynamic contribution because of coupling to 
the vacuum modes. In addition, impurity scattering and 
phonon emission also have contributions to the decoher- 
ence. However, in principle, their effects can be min- 
imized by more precise fabrication technology and by 
cooling the system. 



VII. DYNAMICAL LOCALIZATION OF TWO 
INTERACTING ELECTRONS 

In section IV we have shown that a trivial way of main- 
taining localization is to suddenly shift the constant field 
to another value once the electrons are in the localized 
state \RR). In this section we study the possibility of re- 
maining the localization with the oscillatory electric field 
F(t) = Fismujt Because the localized state \RR) 
may be always produced from the ground state as de- 
scribed above, we therefore suppose in the following dis- 
cussions that the system starts with the localized state 
\RR). 

In the presence of time-dependent electric field, the 
evolution of the system can not be solved in a closed 
form because [Hi(ti), Hi(t 2 )] 7^ 0. However, time pe- 
riodicity of the Hamiltonian (11) enables us to describe 
the dynamics within the Floquet formalism. We numeri- 
cally integrate the motion of equation for time evolution 
operator 



i§- t U(t,0) = #!(*)[/(*, 0), 



(32) 



G 



and diagonalize U(2tt/u),0) to obtain the quasienergies 
{e a ,i \ and Floquet states {| u Q ,i(0))} at time t = 0. Here 
the quasienergies e a ^i are confined to the first Brillouin 
zone and, at F\ = 0, connected to E a + Iuj. The index I 
counts 'how many photons' have to be subtracted from 
the unperturbed energy level E a in order to arrive in the 
first Brillouin zone. The Floquet state | u a ,i(t)) can be 
obtained from the eigenvalue equation 

(Hi(t) - i^) | u a ,,(t)) = e a .i | u a ,j(i)>, (33) 

where a — 1,2,3. Note that the Hamiltonian (11) re- 
mains invariant under the combined spatial reflection and 
time translation t — > t + ir/uj. An immediate consequence 
of this dynamical symmetry is that each Floquet state is 
either odd or even p8| , ^9| . When the driving amplitude 
is switched off adiabatically, F\ — > 0, the Floquet states 
are connected with the stationary eigenstates in Eq. (15) 
as follows [|lj 

I u a .i(t)} ->| u° a /t)) = <p a exp(ilLut). (34) 

Thus we can easily determine the dynamical parity of the 
Floquet state | u a< i(t)). 

We present in Figs. 9(a)-(b) the quasienergies ver- 
sus the amplitude F\, where the values of driving fre- 
quency uj are respectively chosen to satisfy u = 3.1a; 
and u — 2ui, respectively. In Fig. 9(a) we see that the 
quasienergies £2,-3 and £3.-3 with different parity form 
an exact crossing at Fi — 1.05kV/cm. When F\ — > 0, 
all three quasienergies arrive at the unperturbed cases 
by e% = E u e%_ 2 = E 2 - 3lu, £§, _ 2 = E 3 - 3co. For 
21oj — > Ei — Ei, £i i0 and £3,-2;, belonging to different 
parity, are also allowed to cross, as shown in Fig. 9(b) 
where the value of index is I — 1. In the exact same way, 
when (21 + 1)lu — ► E 3 — E\ two quasienergies £1,0 and 
£3 _( 2 ;+i) can develop into the crossing for special value 
of the amplitude F\ (not shown here). 

To elucidate the effect of the exact crossing on the 
quantum mechanic behavior of the system we present in 
Fig. 10(a) the time evolution of PRRit) with the sys- 
tem parameters corresponding to the first crossing be- 
tween £2,-3 and £3,-3 shown in Fig. 9(a). For compar- 
ison we also plot Prr(£) in Fig. 10(b) for the value of 
amplitude F\ — 0.5kV/cm. In Fig. 10(a) we can see 
that during time development the probability Prr(1) re- 
mains near 1 as if the two electrons were frozen in the 
right dot. Thus at the exact crossing of £2,-3 and £3,-3 
the dynamical localization builds up although the strong 
Coulomb repulsion between the two electrons prevents 
the system behaving so. Note that the value of u used 
in Fig. 10(a) corresponds to the unperturbed energies 
Ei = -0.016meV, E 2 = 5.6meV, and E 3 = 5.616meV. 
Because E\ is much lower than E 2 and E 3 , the un- 
perturbed eigenstates (p 2 and tp§, both of which have 
a very small component of the delocalized two-particle 
state I LR), become comparable to a doublet in a sym- 
metrical double-dot system. So it is expected that at the 



crossing of the quasienergies £2,-3 and £3,-3 the initial 
localized state can be approximated by a superposition 
of degenerate Floquet states | 1/2,-3(0)) and | u 3i _ 3 (0)), 
which remains localized in perpetuity. This looks like 
the case of a single-electron, two-level system consisting 
of the lowest symmetric and antisymmetric states of the 
double-trap potential, where perfect localization can be 
achieved at the exact crossing between the two Floquet 
quasienergies. It is numerically found that even if the 
Coulomb interaction u is strong enough, the dynamical 
localization can still occur so long as the quasienergies 
£2,™ and e 3)Tn cross each other. Moreover, the value of 
the system parameter 2(ii/u> corresponding to the first 
crossing is about 2.4, which is a root of the zero-order 
Bessel function, suggesting that in this situation the dy- 
namical localization can be approximated by the driven 
two-level model. If the system parameters deviate from 
the level crossing, then the dynamical localization ceases 
to exist and Prr oscillates between and 1 in the time 
development, as shown in Fig. 10(b). 

More surprisingly, we notice that compared with non- 
interacting (u = 0) case, the presence of Coulomb repul- 
sion enhances the degree of localization to a larger extent. 
To illustrate this feature we present in Fig. 10(c) time 
evolution of Prr with the parameters corresponding to 
the first exact crossing of the quasienergies in the absence 
of Coulomb interaction. It shows in Fig. 10(c) that al- 
though dynamical localization maintains during time de- 
velopment, the degree of localization is much lower than 
that shown in Fig. 10(a), suggesting Coulomb-enhanced 
localization. From Eq. (15) we know that due to strong 
Coulomb interaction, the initial localized state \RR) can 
be approximated by a superposition of two eigenstates 
W { 2 A) ) and |v4 S) > as \RR) ~ (l/y/2){\ V f\ + |^)). 
Thus from what we have learned in single-particle double- 
trap system, it is not difficult to understand why the 
localization can be dynamically remained at the exact 
crossing of two Floquet states developed from \<p 2 A ^) and 

|</?3 ). Similarly, in the absence of Coulomb interac- 
tion, the initial localized state \RR) can be approximated 
by a superposition of all three eigenstates as \RR) = 

(l/V2)(\ip[ S) ) + \(p 2 )) + \<p ( 3 S) ), suggesting dynamical 
localization at the crossing among three quasienergies. 
However, the fundamental difference lies in the fact that 
when the strong interaction is present, the tunneling cou- 
pling is (LL\H\\RR) = E3 — E 2 ~ Aw 2 /u. In the absence 
of Coulomb interaction, whereas, the tunneling coupling 
is (LL\Hi\RR) = (1/2) (E 1 + E 3 ) - E 2 = 2w. Thus we 
can see that the coupling between the two localized states 
greatly decreases in the presence of strong Coulomb in- 
teraction, which leads to fundamental increase of local- 
ization degree. 

We turn to study the dynamics of the system at the 
crossing of the quasienergies £1,0 and £2,-2 m Fig. 9(b), 
using the same initial state condition. The result is shown 
in Fig. 10(d). In contrary to that shown in Fig. 10(a), 
it shows in Fig. 10(d) that at the level crossing of £1,0 
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and £2,-2 dynamical localization does not happen and 
Prr (f ) oscillates between and 1. Note that the level 
crossing of £x,o and £2.-2 induces strong participation of 
the Floquet state | Ui t o(t)} during the time evolution of 
the system, and the most component in | uifi(t)) is the 
delocalized two-particle state | LR) . Therefore the strong 
mixture of the Floquet state | 1x1,0 (t)} in the evolution 
of the system will lead to complete destruction of the 
dynamical localization, as shown in Fig. 10(d). 

In the above discussions we have ignored higher-lying 
single-particle states; this requires that the frequency of 
the external field is much lower than single-particle level 
spacing. In the presence of decoherence due to environ- 
mental dissipation, long dephasing time should be re- 
quired. A detailed analysis of the effect of a decohering 
environment will be given elsewhere. 



VIII. CONCLUSIONS 

In conclusion we have shown how dynamical localiza- 
tion and entanglement of two interacting electrons in a 
double quantum dot system can be accessed by external 
electric fields. We have found that (i) The presence of 
a constant electric field may induce the complete Rabi 
oscillation between the delocalized state | LR) and lo- 
calized state I RR). Thus Starting from the delocalized 
ground state, we can prepare a fully localized state. The 
localization can be maintained by switching the field to 
another nonzero value, (ii) The two electrons oscillate 
between the delocalized state and two localized states in 
the presence of a resonant oscillatory field. With the os- 
cillatory field turned off at time when the probabilities 
for finding the electrons in the left and right dot are iden- 
tically 0.5, the two electrons remain maximally entangled 
in the subsequent time evolution. Thus a selective pulse 
of oscillatory field can be used to implement maximally 
entangled Bell states in a two-electron two-dot system, 
(iii) Although the Coulomb repulsion is very strong, the 
two initially localized electrons can stay localized dur- 
ing time evolution. It is also shown that compared to 
non-interacting case, the Coulomb interaction enhances 
the degree of localization to a larger extent. We expect 
the present results are useful in exploiting the coherent 
control of electrons in quantum dot systems. 

ACKNOWLEDGMENTS 

This work was supported partly by the National Nat- 
ural Science Foundations of China under Grant No. . 



[3] R. Bavli and H. Metiu, Phys. Rev. Lett. 69, 1986 (1992). 
[4] J. M. Gomez Llorente and J. Plata, Phys. Rev. A 45, 
R6958 (1992). 

[5] F. Grossmann and P. Hanggi, Europhys. Lett. 18, 571 
(1992). 

[6] M. Holthous and D. Hone, Phys. Rev. B 47, 6499 (1993). 
[7] R.I. Cukier and M. Morillo, Chem. Phys. 183, 375 (1994). 
[8] T. Zuo, S. Chekowski, and A.D. Bandrauk, Phys. Rev. A 

49, 3943 (1994). 
[9] H. Pu, S. Raghavan, and N. P. Bigelow, Phys. Rev. A 61, 

023602 (2000). 

[10] P. Zhang, A.Z. Zhang, S.Q. Duan, and X.-G. Zhao, Phys. 

Lett. A 287, 105 (2001). 
[11] P. Zhang, Q.-K. Xue, X.-G. Zhao, and X.C. Xie, Phys. 

Rev. A 66, 022117 (2002). 
[12] P.I. Tamborenea and H. Metiu, Europhys. Lett. 53, 776 

(2001). 

[13] A. Steane, Rep. Prog. Phys. 61, 117 (1998). 

[14] Q.A. Turchette, C.S. Wood, B.E. King, C.J. Myatt, D. 

Leibfried, W.M. Itano, C. Monroe, and D.J. Wineland, 

Phys. Rev. Lett. 81, 3631 (1998). 
[15] E. Hagley, X. Maitre, G. Nogues, C. Wunderlich, M. Brune, 

J.-M. Raimond, and S. Haroche, Phys. Rev. Lett. 79, 1 

(1997) . 

[16] A. Rauschenbeutel, G. Nogues, S. Osnaghi, P. Bertet, M. 
Brune, J.-M. Raimond, and S. Haroche, Science 288, 2024 
(2000). 

[17] CA. Sackett, D. Kielpinski, B.E. King, C. Langer, V. 
Meyer, C.J. Myatt, M. Rowe, Q.A. Turchette, W.M. Itano, 
D.J. Wineland, and C. Monroe, Nature 404, 256 (2000). 

[18] F. R. Waugh, M. J. Berry, D. J. Mar, R. M. Westvelt, K. 
L. Campman, and A. C. Gossard, Phys. Rev. Lett. 75, 705 
(1995); C. Livermore, C. H. Crouch, R. M. Westvelt, K. L. 
Campman, and A. C. Gossard, Science 274, 1332 (1996); 
R. H. Blick, D. Pfannkuche, R. J. Haug, K. v. Klitzing, and 
K. Eberl, Phys. Rev. Lett. 80, 4032 (1998); R. H. Blick, 
D. W. vander Weide, R. J. Haug, and K. Eberl, ibid. 81, 
689 (1998). See also references in these publications. 

[19] For recent articles on double quantum dots with few elec- 
trons see: Y. Tokura, D.G Austing, and S. Tarucha, J. 
Phys.: Condens. Matter 11, 6023 (1999); B. Partoens, 
A. Matulis, and F.M. Peeters, Phys. Rev. B 59, 1617 
(1999); D.G. Austing, T. Honda, K. Muraki, Y. Tokura, 
and S. Tarucha, Physica B 249, 206 (1998); and references 
therein. 

[20] B.E. Kane, Nature 393, 133 (1998). 

[21] D. Loss, and D.P. DiVincenzo, Phys. Rev. A 57, 120 

(1998) . 

[22] P. Zhang and X.-G. Zhao, Phys. Lett. A 271, 419 (2000). 

[23] A. Imamoglu, D.D. Awschalom, G. Burkard, D.P. DiVin- 
cenzo, D. Loss, M. Sherwin, and A. Small, Phys. Rev. Lett. 
83, 4204 (1999). 

[24] L. Quiroga and N.F. Johnson, Phys. Rev. Lett. 83, 2270 

(1999) . 

[25] L. Jacak, P. Hawrylak, and A. Wojs, Quantum dots 

(Springer, Berlin, 1997) 
[26] G. Burcard, G. Seelig 



[1] M. Grifoni and P. Hanggi, Phys. Rep. 304, 229 (1998). 

[2] F. Grossmann, T. Dittrich, P. Jung, and P. Hanggi, Phys. 
Rev. Lett. 67, 516 (1991). 



mat/9910105) 



and D. Loss, preprint (cond 



[27] J. Schlicmann, D. Loss, and A.H. MacDonald, preprint 
(|cond-mat/000908^). 



8 



[28] A. Peres, Phys. Rev. Lett. 67, 158 (1991). 

[29] J. Y. Shin and H. W. Lee, Phys. Rev. E 53, 3096 (1996). 



corresponding to the exact level crossing shown in Fig. 
9(b). 



Figure captions 

Fig. 1 Sketch of two interacting elctrons in a coupled 
quantum dot driven by electric fields. 

Fig. 2. Time evolution of the probabilities Plr (solid 
line), Pll (dashed line), and P RR (dotted line) for the 
value of the strength of constant electric field satisfying 
resonant condition no = u. 

Fig. 3. (a) Electric field that imposes on the double 
quantum dot system; (b) Time evolution of the proba- 
bilities Plr (solid line) and Prr (dotted line) under the 
influence of the electric field shown in (a). The other 
system parameters are the same as that used in Fig. 2. 

Fig. 4. (a) Electron number distribution of right 
dot in the ground state as a function of the strength of 
constant electic field; (b) Energy spectrum of the driven 
two-electron system as a function of the strength of a 
constant electric field. 

Fig. 5. (a) Electric field that imposes on the double 
quantum dot system; (b) Time evolution of the proba- 
bilities Plr (solid line) and Prr (dotted line) under the 
influence of the electric field shown in (a); (c) Time evo- 
lution of the probability pBeii for finding the maximally 
entangled Bell state {<p = — ?r/2) in the presence of the 
Coulomb interaction (solid line) and in the absence of the 
Coulomb interaction (dotted line). 

Fig. 6. Quasienergy spectrum of the driven two- 
electron system as a function of the frequency of oscilla- 
tory electric field. 

Fig. 7. (a) Time evolution of the probabilities Plr 
(solid line), Pll (dashed line), and Prr (dotted line) 
for the value of the frequency of oscillatory electric field 
lo = u; (b) Time evolution of the probability pBeii for 
finding the maximally entangled Bell state (<f) — tt) , The 
system parameters are the same as that used in (a). 

Fig. 8. (a) Electric field that imposes on the double 
quantum dot system; (b) Time evolution of the probabili- 
ties Plr (solid line), Pll (dashed line), and Prr (dotted 
line) under the influence of the electric field shown in (a); 
(c) Exact numerical (solid line) and approximate analytic 
(dotted) results of time evolution of the probability pBeii 
for finding the maximally entangled Bell state (<f> = tt) 
in a coupled quantum dot. The other system parameters 
are the same as that used in Fig. 7. 

Fig. 9. Floquet spectrum of the driven two-electron 
system as a function of the strength of oscillatory electric 
field for the value of the effective Coulomb interaction (a) 
u = 3.1a;; (b) u = 2uo. 

Fig. 10. Time evolution of the probability P RR (t) 
for four different kinds system parameter values (a) u — 
3.1a; and F\ = 1.05kV/cm, corresponding to the exact 
level crossing shown in Fig. 9(a); (b) u = 3.1a; and 
Fi = 0.5kV/cm, for comparison with the case shown in 
(a); (c) u = 0, u> = 1.81meV, and F\ = 2kV/cm, cor- 
responding to the exact level crossing in the absence of 
Coulomb interaction; (d) u = 2u and F\ = 0.52kV/cm, 
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